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Abstract 



We study nonsingular cosmological scenarios in a general D-dimcnsional string effective 
action of the dilaton-modulus-axion system in the presence of the matter source. In the 
standard dilatonic Brans-Dicke parameter (to = —1) with radiation, we analytically obtain 
singularity-free bouncing solutions where the universe starts out in a state with a finite 
curvature and evolves toward the weakly coupled regime. We apply this analytic method 
to the string-gas cosmology including the massive state in addition to the leading massless 
state (radiation), with and without the axion. We numerically find bouncing solutions which 
asymptotically approach an almost radiation-dominant universe with a decreasing curvature 
irrespective of the presence of the axion, implying that inclusion of the matter source is 
crucial for the existence of such solutions for oj — —1. In the theories with oj ^ —1, it is 
possible to obtain complete regular bouncing solutions with a finite dilaton and curvature 
in both past and future asymptotics for the general dimension, D. We also discuss the 
case where dilatonic higher-order corrections are involved to the tree-level effective action 
and demonstrate that the presence of axion/modulus fields and the matter source does not 
significantly affect the dynamics of the dilaton-driven inflation and the subsequent graceful 
exit. 



1 Introduction 



String cosmology has continuously received much attention T with the development of string 
theory as a possible candidate to unify all fundamental interactions in nature For cos- 
mologists it is very important to test the viability of string theory by extracting cosmological 
implications from it. In particular string cosmology has an exciting possibility to resolve the 
big-bang singularity which plagues in General Relativity. 

Among the string-inspired cosmological scenarios proposed so far, the pre-big-bang (PBB) 
model ]3|i (3| based on the low energy, tree-level string effective action has been widely stud- 
ied (see ref. [HI for a recent review). In this scenario there exist two disconnected branches, 
one of which corresponds to the dilaton-driven inflationary stage and another of which is the 
Friedmann branch with decreasing curvature. Then string corrections to the effective action 
can be important around the high curvature regime where the branch change occurs. Recently 
proposed Ekpyrotic/Cyclic cosmologies |S1EI have a similarity to the PBB scenario in the sense 
that the description in terms of the tree-level effective action breaks down around the collision 
of two branes in a five dimensional bulk. 

When the universe evolves toward the strongly coupled, high curvature regime with growing 
dilaton, it is inevitable to implement higher-order corrections to the tree-level action. Indeed it 
was found that two branches can be smoothly joined by taking into account the dilatonic higher- 
order derivative and loop corrections in the context of PBB and Ekpyrotic jlOj scenarios. In 
the system where a modulus field is dynamically important rather than the dilaton, Antoniadis 
et al. showed that the big-bang singularity can be avoided by including the Gauss-Bonnet 
curvature invariant coupled to the modulus 

In contrast we can consider the case where the dilaton evolves toward the weakly coupled 
region. In the brane-gas cosmology ^1 El El E| and the recently proposed string-gas 
cosmology |17| . the universe starts out from a dense and hot state filled with a gas of fun- 
damental branes and strings in thermal equilibrium. These are originated from the work of 
Brandenberger and Vafa ^Sl (see also ref. ^HI) in an attempt to construct singularity- free cos- 
mological scenarios using T-duality. In these models it is typically assumed that the cosmic 
evolution does not reach the strongly coupled regime in order to avoid the breakdown of the 
effective ten-dimensional background description (see ref. [SUj for the parameter regions where 
the effective ten-dimensional description of the physics is no longer valid). In particular, in the 
context of the string-gas cosmology, it was shown that only the large 3 dimensions expand and 
asymptotically approach the radiation dominant universe while the small 6 or 7 dimensions are 
kept undetectably small jTT]. 

In order to obtain viable nonsingular cosmological scenarios in string theories, the string 
coupling and the scalar curvature are required to be finite during the cosmological evolution. In 
the standard tree-level action with the dilatonic Brans-Dicke parameter w = — 1, it is generally 
difficult to construct nonsingular solutions in the absence of the matter source term. This comes 
from the fact that the evolution toward the strongly coupled regime reaches the curvature 
singularity in the tree-level action and that the evolution toward the weakly coupled regime 
does not have enough power to lead to bouncing solutions. Although bouncing solutions were 
found in ref. \21\ \'2'2\ l2,Sj by including the axion field in the tree-level action, this corresponds 
to the case where the string coupling and the scalar curvature diverge in both past and future 
asymptotics. Therefore this can not be regarded as ideal singularity-free bouncing solutions 
unless some higher-order effects alter the dynamics of the system when the dilaton enters the 
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strong coupling regime. 

We can take into account the radiation-like matter source which comes, for example, from 
the ideal string-gas or the 5-form of the Ramond-Ramond sector of the type IIB superstring 
|24| I25| I26j . In the presence of the radiation matter and the axionic field in d = 3 spatial 
dimensions, Constantinidis et al. [21] showed that bouncing solutions can be obtained when 
the dilatonic Brans-Dicke parameter uo is negative (i.e., including the co = —1 case). When 
—4/3 < a; < these solutions exhibit divergent behavior of the scalar curvature as the dilaton 
evolves toward the large coupling regime {(p — > oo), in which case Constantinidis et al. did 
not regard them as ideal nonsingular bouncing solutions. Meanwhile these solutions, in another 
asymptotic, tend to approach the radiation dominant universe with a decreasing curvature and 
a finite dilaton. Then it should be possible to have regular bouncing solutions if the initial 
state of the universe is in the weakly coupled region with a finite curvature as in the case of 
string/brane-gas cosmology. 

In this paper we shall investigate singularity-free cosmological scenarios based on the effective 
action of the dilaton-modulus-axion system in the presence of the matter source term. We keep 
the dilatonic Brans-Dicke parameter co arbitrary so that the effective action includes a wide 
variety of theories such as the F-theory [27\ I28j or the multidimensional theory ^2] • Indeed it 
was recently found that in the presence of the axion and the pure radiation complete regular 
bouncing solutions were found in both past and future asymptotics when u; < — 3/2 in the case of 
d = 3 spatial dimensions |25| I26j . With an application to string/brane-gas cosmology in mind, 
we shall extend the analysis to the general d spatial dimensions with/without the axion and 
modulus fields. In fact there exist a number of bouncing solutions even for the tree-level action 
in the presence of the radiation source term. In particular we will make numerical analysis in 
the context of string-gas cosmology with and without the axion by including the massive state 
together with the leading massless state (radiation) and demonstrate the existence of nonsingular 
bouncing solutions even for lo = —1. 

When the system evolves toward the strongly coupled regime, it is likely that the evolution 
is altered by the higher-order corrections to the tree- level action. Therefore we shall also make 
numerical simulations for the case where the dilatonic high-order derivative and loop corrections 
are included in the dilaton-modulus-axion system with a matter source term. We will show 
that these higher-order corrections typically dominate the dynamics of the system, by which the 
big-bang singularity can be avoided as in the case where the only dilaton is present. 



2 Model 

Consider the following effective action with three scalar fields (p, x and a: 



S = d^XyJ-QD 



(1) 



where R is the scalar curvature in Z)-dimensions with metric g^^ with = Aetg^y. Cm 
represents the matter source term and Cc corresponds to the higher-order corrections to the 
tree-level action. The action contains a variety of models derived from string theories. In 
what follows we shall briefly show such examples. 



Low energy effective string action 
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The lowest order D-dimensional string effective action coming from the beta function of 
the string world-sheet is characterized by ^ 



where ip is the dilaton, H^^x is the antisymmetric tensor field. Let us compactify the 
k = D — A dimensions by introducing a modulus field, e^^. We define an effective dilaton 
in four dimensions as 99 = ■0 ~ ^X- We also introduce a pseudo-scalar axion field a as 
fjij.i^>' — ^ip^fiuXp-^ ^a,king into account the fact that iJ^'^-^ has only one degree of 
freedom in 4 dimensions 1 . Then we get the action with D = A, uj = —1, x = V^X) 
n = —1, Cm = and Cc = 0. 

String /Brane- gas cosmology 

The string-gas cosmology |17| is an attempt to explain why and how the three spatial 
dimensions become large and other 6 or 7 are kept to be small. The bulk action in 10 
or 11 dimensions is assumed to take the same form as ((21), but we have additional source 
matters due to the presence of the ideal string-gas in thermal equilibrium. The simplest 
version of the string-gas cosmology [T7j corresponds to w = —1, x = 0; = 0) 7^ 
and £c = in (P). Instead of introducing a modulus field, one considers the evolution 
of one scale factor for [D — 1) spatial dimensions or two scale factors for large and small 
dimensions jl7j . In this work we shall also analyze the case where the axion is present 
(cT 7^ 0). Note that the brane-gas cosmology discussed in refs. ^lEI belongs to this class 
with a different matter source when we use the low energy bulk action with a; = — 1. 

Cosmology with p-brane or Dp-brane solitons 

The p-brane or Dp-brane are the solitonic degrees of freedom on which string endpoints 
live. These can be fundamentally important when the string coupling is large, since they 
are light in the high curvature regime j2Ul l3Uj . In the string u-model whose metric is 
minimally coupled to the p-brane, Duff et al. [HJ showed that the effective action can be 
described by ^ with x = 0, a = 0, Cm 7^ and 

{D-l){p-2)-p^ 

(Z)-2)(p-2) -p2 ' ^""^ 

where p = p + 1. From this we have to = —1 for the string (p = 1). However lo depends 
upon the values of d and p for the p-branes with p 7^ 1. Cosmology in this scenario was 
investigated in refs. |321 133j in the absence of higher-order corrections to the tree- level 
action (see also ref. [Hi])- 

F- theory 

The superstring type IIB theory can be reformulated geometrically within the framework 
of 12 dimensional theory- called i^'-theory [271 [28) . After the dimensional reduction from 
12 to 10 dimensions, the following action may be obtained |26j : 



^=2 



i J d'^x./^ [e-'^ (^i?io + 3(Vv9)2 - -Lf^^.F'^-^) - ^F^^F'^-' 



(4) 
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where the terms F^^ and F^^^x come from the Ramond-Ramond (RR) sector of the type IIB 
superstring. This action is different from the standard string theory with lo = —1. Making 
the similar dimensional reduction from 10 to 4 dimensions as @ with the isotropization 
of the RR terms, we get the action with uj = —3, x 7^ Oi <7 7^ (n = — 1) and the 
radiation fluid term {Cm 7^ 0) coming from the 5-form of the RR sector. The axionic 
terms minimally coupled to the dilaton can also appear from the same RR sector |26) . 
This corresponds to the case of adding another field with n = in the action 

• Multidimensional cosmologies 

The low energy limit of the string effective action may be reformulated in the context of 
multidimensional theories j^Hj . In its simplest form the multidimensional theories contain 
only the geometry, in which case one has uj = {l — d)/d with x = cr = by the dimensional 
reduction (d is the number of compactified dimensions). When a 2-form field or a conformal 
gauge field with (d + 4)/2-form is present, we have an axionic term a in eq. with 
n = —2/d + 1 or n = —2/d, respectively. Therefore we keep the value of n arbitrary so 
that multidimensional theories are also included in (^. 

Hereafter we shall analyze the cosmological solutions based on the action . For convenience 
we call the fields if, x and a as dilaton, modulus and axion, respectively, unless otherwise 
specified. We do not include the potential of scalar fields in the action^. Extending the analysis 
of refs. j24| I25| I26j with D = 4, we will study the cosmological evolution in the general D- 
dimensional action with a time-dependent modulus field (x 7^ 0). In applying to string-gas 
or p-brane cosmologies, it is convenient to consider the case of the general dimension D. In 
addition the presence of the internal space (modulus) can alter the parameter range where the 
nonsingular solutions exist. We should also keep in mind that higher-order corrections to the 
tree-level action are inevitably important in the strongly coupled regime (e'^ ^1)- later 
sections we shall numerically investigate the dynamics of the system in the action when the 
dilatonic a' and loop corrections are taken into account (£c 7^ 0). 

Let us consider the case where the D = d + 1 dimensional spacetime is described by a flat 
Priedmann-Robertson- Walker (FRW) metric with the line element 



where a{t) is the scale factor. In the context of the string (brane)-gas cosmology, it is possible 
to consider the case with two scale factors corresponding to the "large" and "small" dimensions. 
It was shown in ref. [T7] that only the large dimensions can be dynamically important whereas 
the small dimensions are kept to be nearly constant. In such cases it is sufficient to analyze the 
case of the one scale factor described by eq. ©. However we will also discuss the case of two 
scale factors in sec. 5. 

Defining a new scalar field. 



ds 



,2 



dt'^ + a^{t)6ijdx^dx^ , 



(5) 



(6) 



the variation of the action yields the following equations of motion 




See ref. |35| for the construction of nonsingular cosmological solutions including the dilaton potential. 
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H = -jj^Yj^\^i^^ + (l)X^ + (l)'"&^-H<j) + 4> + P + P + Pc+Pcj , (8) 

4> + dH<p - ^ (u^ +x^ + ncP^'-^a^ + i? + A^^ = , (9) 

X+(^dH + ^^X = 0, (10) 

a+ (^dH + a = 0, (11) 

p + dH{p + p) = 0, (12) 

where H = d/a is the Hubble rate. Here p and p stand for the energy and pressure density of the 
matter source, respectively, pc and pc correspond to the dilatonic higher-order corrections to the 
tree-level action with stress-energy tensor T/f = {—pc^Pc^Pc^Vcj- comes from the variation of 
Cc with respect to </>. 

The scalar curvature, i?, is expressed as 

R = d{d+l)H'^ + 2dH (13) 

= -co^-X^-r-'^' + j^^Y^{-d'H^-d4> + p-dp + p,-dp,). (14) 

Substituting this relation for eq. Q, we get 

I , ^^i , {l-n){d-l)r(T^-2{p-dp)-2{p,-dp,)-{d-l)cl>A, 

^ + 2Md-i) + 4 = 

The equations (fTUIl and (fTT|) are easily integrated to give 

where A and B are integration constants. Making use this relation and introducing a new time 
parameter, 

-d 



T = J a-'^dt , (17) 

the (p equation (|15j) can be written in the form: 

d^ (l_n)(d-l)S2 _ 2(p-(ip) + 2(pe-rfpc) + (d-l)0Ae 2d ^ . ..ox 
dr2 ^ 2V((i - 1) + 4 2V(d-l)+4 ^ ^ 

In the absence of the correction Cc, the last term of eq. (|18|) vanishes for p = p = (no matter 
source) or p = dp (radiation). In these cases it is possible to obtain analytic solutions for the 
general spatial dimension d even in the presence of axion and modulus fields. Notice that the 
simplest pre-big-bang scenario [S] corresponds to the case of p = p = without the axion field. 
The radiation matter (p = dp) can appear naturally in the context of string-gas cosmology 
at finite temperature ^7] or the 5-form existing in the RR sector of the type IIB superstring 
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When p = dp, the integration of eq. H12() gives p = poa ('^+^), with pQ being a constant. 
Introducing a new parameter, b = a(f?'^^'^~'^\ eq. can be written for £c = as 

d{d - l)4>^ ' = 2^06^^-^ +A' + B'^'-- + ^id-l)+d ^„ ^ ^^g^ 

where a prime denotes the derivative with respect to r. This equation includes the case without 
the matter source by setting po = 0. 

In the next two sections we shall derive the analytic solutions for the scale factor and the 
dilaton in the case of the radiation or no matter source without the correction Cc- For the general 
equations of state with p = jp, it is not so easy to find analytic solutions due to the fact that 
the last term in eq. (|18)) does not vanish^. In sec. 5 we will analyze the dynamics of string-gas 
cosmology as an application by including the massive state (Kaluza-Klein and winding modes) 
in addition to the leading massless state (radiation), with and without the axion. In this case 
the string-gas state is not purely the radiation. Since it is difficult to proceed analytic approach 
when the higher-order correction Cc is present, we shall make numerical analysis separately in 
sec. 6 for the general action Hereafter we will consider the case of d > 1. 



3 Low energy tree-level action without axion (cr = and Cc = 0) 



Let us first analyze the low energy tree- level action without the axion (Cc = and a = 0). 
When p = p = (no matter source) or p = dp (radiation) , eq. ((TH|l is integrated to give 



cr . 



(20) 



where c is a constant. We can set another integration constant to be zero without loss of 
generality. Since (j) = e"*^ is required to be larger than zero, we have r > for c > and r < 



for c < 0. The string coupling, gj 



diverges as r — > 0. 



Since we are now considering the case with B = 0, eq. (fT!l|) yields 



d6 



where 



M 



2/90 



d{d-l) ' 



N 



d{d - 1) 



dr 



■^2 ^ i^{d-l) + d ^^ 



d-1 



(21) 



(22) 



The integral in the l.h.s. of eq. ()21() is different depending on the sign of A'^. For example, N is 
positive in the standard string theory (oj = —1) with d > 1. However, can be negative for the 
theories such as the solitonic p-brane/Dp-brane or multidimensional theories. In what follows 
we shall discuss the positive and negative cases separately. 

^Note, however, that it is stiU possible to obtain analytic solutions for some special equations of state in the 
absence of the axion and the modulus 1321. 
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3.1 Case of > 



When > the integration of eq. (|21|) leads to 

6^-1 1 



(23) 



where 



i = io{cTf With q = ±^{d-l). (24) 

2c 

Here is a positive integration constant. Then the analytic solution for the scale factor can be 
written as 

= /(<^)'/("-') with /(</)) = , '^-^ —. (26) 

eo<A^+^ - (M/Co)^^-'^ 

This indicates that the evolution of the scale factor is determined in terms of the field (f) for the 
case (jSnj- The simplest pre-big-bang scenario with a constant modulus corresponds to the case 
of ^ = 0, I? = 0, M = 0, (i = 3 and oj = —1, thereby yielding q = ±^^3/6 and 

a = ^^(-3±V3)/6 (27) 

This represents an inflationary solution with decreasing (j) from the weakly coupled regime (</> 
oo) to the strongly coupled regime {(f) 0). Note that the solution eventually reaches the 
curvature singularity as — > 0. If we consider the case of increasing (j), the solution presented 
in eq. ()27() describes a contracting universe. In this case the solution does not connect to our 
expanding branch. 

The situation changes when the radiation is taken into account (M ^ 0). Differentiating 
eq. with respect to (j), we have 

d/.,. , ^ (1 + 2q)Coc^'2+i - (1 - 2g)(M/Co)0^'? 

(^0^^ - (M/Co)^-)^ • ^''^ 

This means that it is possible to have bouncing solutions with the increase of (p as long as 
— 1/2 < q < 0. In this case, if the condition (1 + 2g)^o</'2+'? > (1 - 2q){M/(,o)<P^''' is satisfied 
initially with a finite value of the dilaton ((pi ^ 0), the universe contracts at the initial stage. 
This is followed by a bounce due to the growth of the (1 — 2q){M/^Q)<f>2~'3 term. The bounce 
occurs when the field (j) passes at 

The universe begins to expand for with the increase of (j). As found from eq. ()26|) the 

scale factor goes to infinity as (j) approaches the critical value, 

(30) 



(29) 
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Let us find the asymptotic behavior of a around (p = (f)c- Setting (f> = (j)c — cAr with At being 
a small time parameter, we find the scale factor (|26|) is proportional to (Ar)~^/^'^~^^. Since the 
cosmic time t is related with r by the relation dt/dr = a'^, the evolution of the scale factor is 
described by 

^^^2/(d+l) (f^^, (31) 

which shows the radiation dominant behavior. 

In the limit of ^ — >■ 4>c, we have t — >■ oo as At 0. Then the scalar curvature asymptotically 
approaches zero as i? oc for cp ^ cpc- In addition the string coupling, = is also finite 
in this case. This means that nonsingular bouncing solutions with finite curvature and dilaton 
can be obtained by taking into account the radiation in the low energy tree-level string action. 
Notice that we are not considering the regime around the strong coupling limit ((/>—> 0). The 
universe is assumed to start out from the region with finite dilaton and curvature, as in the 
string/brane-gas cosmologies. Then the dilaton evolves toward the weakly coupled region (see 
fig.^, which asymptotically approaches the radiation dominant universe given by ^'M\i . In sec. 5 
we will numerically investigate the existence of such bouncing solutions in string-gas cosmology 
taking into account the massive state in addition to the leading massless state. 

These bouncing solutions exist for — 1/2 < g < 0, whose condition is written as 

oj<-^. (32) 

In the absence of the modulus, this condition yields c<j < 0. Note, however, that the large initial 
value of the modulus restricts the range for the existence of the bouncing solutions. 

If the condition, (l+2g)^o</'^"^^ < (l-2g)(M/^o)</'^"^, is satisfied initially with -1/2 < g < 0, 
the scale factor (|26|) grows from the beginning with the increase of (j). Then it asymptotically 
approaches the radiation dominant solution given by eq. H31|) . This solution was numerically 
found as well in the context of string- gas cosmology [l^- Note that when \q\ > 1/2 we have 
either ever expanding or contracting solutions. 

It is possible to obtain bouncing solutions provided the condition H32|) is satisfied and the 
dilaton decreases. In contrast, it is worth mentioning the case where the dilaton evolves toward 
the strongly coupled regime. In this case, the scalar curvature does not necessarily converge in 

1+29 

the limit of — > 0. Around the region of = 0, the scale factor is estimated as a oc (/> ^-i for 
— l/2<g<0. Combining this relation with dt/dT = a'^, we have 

l + 2dq 

t = to + CT , (33) 

where to and c are constants. Then the scale factor around = is expressed as 

l + 2q 

a OC (t - to) . (34) 

When q > —l/{2d), one has t — > oo as r ^ from eq. Since the power in eq. (|5H) is 

positive in this case, the universe expands for r ^ 0. Meanwhile we have t — > to r — > for 
q < —l/{2d). Since we are considering the case of g > —1/2, the power in eq. ()34() is negative 
for —1/2 < q < — l/(2d), thereby leading to the growth of a as t ^ to- 

Since the scalar curvature evolves as i? oc (t — to)~^ for the solution (|34j) . it asymptotically 
decreases toward zero only for q > —l/{2d). When —1/2 < q < —l/{2d), R diverges in the 
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Figure 1: A nonsingular bouncing solution in the low energy tree- level action without the axion 
in the presence of the radiation (p = dp). With the growth of (/> = e~'^, the universe contracts 
and exhibits the bounce at (j) = (p^,. After that the scale factor begins to grow and the system 
asymptotically approaches the radiation dominant universe as (j) ^ (pc- 



limit of (/> — > 0. Therefore, when the system approaches the strongly coupled regime, we have 
bouncing solutions with finite curvature for q > —l/{2d), i.e., 



1 

d 



^ < _1 _____ . (35) 



This does not include the string case {uj = —1). Combining ()35() with the condition > 0, we 
have — — ^ < uj < — 1 — ^ — When d = 3 and A = (no modulus) , this condition yields 
-3/2 < u; < -4/3. 

We have to keep in mind that the description of the tree-level action ceases to be valid as 
(p evolves toward the strong coupling region {4> ^ 0). The higher-order corrections to the tree- 
level action are expected to be important in this regime. On the contrary, nonsingular bouncing 
solutions toward the weakly coupled regime shown in fig. ^ can be regarded as ideal ones. In 
addition the dilaton and the scalar curvature do not diverge, provided that the universe starts 
out from the region with finite (p and R. 

3.2 Case of < 

Let us next consider the case of A < 0. For example, when d = 3 and A = 0, this condition 
corresponds to a; < —3/2. The only difference relative to the A > case is the integral in the 
l.h.s. of eq. (|21|). For negative A we have 



arctanW — =\og^Q\cT + 4)q\'' with g = ± — - — (d — 1) , (36) 
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where is a positive constant. Then the scale factor can be written as 



\N\ 1 

M V^cos(logCo0«) 



2/{d-l) 



(37) 



This indicates that nonsingular bouncing solutions exist in the interval 



1 

— exp 



< 



< -p- exp 



2. + i 

•' 2 



7T 



(38) 



where j is an integer. From eq. p7j) we find that the bounce occurs at cos(log Co^*) 
i.e., 



Co 



1/9 



exp 



1 

-arc cos 

iq 



q+y/q' + 1, 



(39) 



Then we have q < from the condition, | cos(log^o<^*)| ^ 1- The asymptotic behavior of the 
scale factor is given by a oc 

f2/(d+i) ^j^g dilaton approaches (p ^ (l/^^"^) exp[(2j ± l/2)7r/g]. 
In addition it is easy to show that the scalar curvature asymptotically decreases in the limit of 
^ (1/Cy'^)exp[(2j ± l/2)7r/g]. Therefore the bouncing solution (|37|) is characterized by the 
radiation dominant universe in both asymptotics with finite values of the dilaton and the scalar 
curvature. We have a sequence of solutions (j^Hj) corresponding to the each integer value of j. 
This is different from the case of iV > where singularity-free solutions with finite dilaton exist 
only for (pi < (p < (j)c {4>i is the initial value). The key point is that the dilaton can be controlled 
without entering the strongly coupled region with ~ 0. 

In the absence of the axion with p = dp (or p = p = 0) considered above, the dilaton 
evolves as the free massless field {<j)" = 0). However the situation changes when the axion is not 
dynamically negligible. 



4 Low energy tree-level action with axion {a ^ and Cc = 0) 



In this section we implement the effect of the axion and discuss the difference of the dynamics 
compared to the case of u = 0. First, the evolution of the dilaton is altered by the existence 
of the second term in eq. H18|) . When Cc = and p = dp (or p = p = 0), we find a conserved 
quantity, E, from eq. (|TH|) : 



2^ 



id-l)B^ 



E 
2" 



2{uj{d - 1) + d}^ 

Making use of this relation, eq. is reduced to the integral form (|21j) with 



M 



2po 



d{d - 1) 



N 



1 



d{d-l) 



a2 ^{d 



l)+d 



d-1 



E 



(40) 



(41) 



The only difference relative to the a = case is that in eq. 1)22(1 is replaced for E. 

The integral in the l.h.s. of eq. ()21() can be performed as in the previous section depending 
on the sign of N. In order to integrate the r.h.s. of eq. (|2H) . it is required to make further 
classifications depending on the sign of the second term in eq. (|4U|) . 
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4.1 uj>-d/{d-l) 

When uj > —d/{d — 1), it is convenient to introduce a new time parameter, 77, satisfying 



^ = ^o(sinr?)^/(^-") with <^o= r^T~\l. ] ■ (42) 



fE{ujid-l)^d}^^/^^-''^ 

From eq. (j^O)) we find the relation d^/dr = zbV^cosjy. Together with eq. (|lT|l we obtain 

dr 2 i+n 

;r = -7^</'o(sin??)i-" . (43) 

dr? (l-n)V£^ 

Making use of this relation it is easy to perform the integral (|2H) . Notice that we are now 
considering the case with > since this is automatically satisfied for uj > —d/{d— 1). The 
scale factor is expressed by the form (|25|) with 

?7 

tan- 



e = Co 

Namely the explicit form of a can be written in terms of rj : 
( 2VN 




1 = ^^\H^- (44) 



2/(d-l) 



^)i/(i-n) {^tanr]/2\i - (M/^o)! tan 7?/2|-'2} 



sm 



(45) 



The evolution of the scale factor is well understood by considering the behavior of and for 

< ?7 < TT (or the region 2j7r < r] < {2j + 1)tt with j being an integer). 

Let us first investigate the case of n < 1 and q < 0. It is sufficient to analyze negative values 

of q only, since the function ^ is symmetric with respect to rj = 7r/2 for the change of q ^ —q. 

As found from eq. 1)251) . ^ is larger than \/M. It is clear from eq. ()25|) that the universe expands 

as C approaches = VM (note that (f) is finite in this limit). Consider the asymptotic behavior 

of the scale factor when ^ evolves toward = Co |tan(r7c/2)|'' with M 7^ 0. Setting r] = r]c — At/ 

2 

in eq. ()25)) . we find that the scale factor evolves as a oc (Ary) <*-i around = ??c- Making use of 

the relation ()45() and dr/dt = a"'^, we have o oc with t ^ 00 as Arj 0. This represents 

the radiation dominant universe with decreasing curvature {R 0) as is similar to the case 

In another limit ry — > 0, one has (/> — > and C ^ 00 for n < 1 and q < 0. Since the evolution 
of the scale factor is approximately described by a oc ^^(1^"*"^) around = 0, we get an 
expanding solution for q > —j^- From eq. ()41() and 1)44(1 this condition yields 

^<-— , (46) 

which is similar to eq. (|32() . Following the similar process by which eqs. H33|) and (|34|1 are derived, 
we get the evolution of the scale factor around r/ = as 

a OC (t - to) i+'^'ti-") . (47) 

Here t = to + cr]~'^^'^~^'^'^^ with to and c being integration constants. It is easy to show that 
the scalar curvature decreases toward zero for g > — ^^^^^^^ , whereas R diverges for (7 < — ^^-^^^^ . 
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If we want to construct bouncing solutions where the dilaton evolves toward (f) = with finite 
curvature, this requires the condition q > — ^(^-n) ' 

.<-l-i-|. (48) 

which is analogous to eq. pSf) . Combining this with the condition u; > —d/{d — l), we have 
—3/2 < u; < —4/3 for d = 3 and A = (no modulus), as shown in ref. |24j . Notice that the 
presence of the modulus makes the allowed range of uj narrower. 

Again we should caution that the string coupling, g"^ = diverges in the asymptotic limit 
0^0. In this sense these solutions can not be regarded as ideals ones, unless some higher- 
order effects are taken into account in the strong coupling regime. Meanwhile when the universe 
begins from a state with finite dilaton and curvature satisfying (pi > 1, it is possible to have 
nonsingular bouncing solutions provided that the solutions approach the region around C = Cc- 
In this case the dilaton is finite for the range 4>i < (f) < (pc with = (po {smr]c)'^^^^~^\ The 
system asymptotically approaches the radiation dominant universe with decreasing curvature. 
One different point relative to the a = case is that the axion tends to work for the dilaton to 
turn back toward the strong coupling regime. Nevertheless the string coupling does not diverge 
as long as the dilaton approaches the value (j) = (j)c instead of = 0. 

When n > 1, both ^ and r] diverge either in the limit of ?? — > vr or — > 0, depending on the 
sign of q. Since the asymptotic evolution of the scale factor exhibits the contracting behavior, 
it is not possible to have bouncing solutions for n > 1. 

At the final of this subsection, it is worth mentioning the case without the matter source 
(M = 0). When n < 1 and q < 0, the scale factor approaches the asymptotic solution given 
by eq. 1)471) in both limits of r/ ^ and r/ — > vr. Therefore if the condition (|48() is not satisfied, 
the scalar curvature diverges in both past and future asymptotics together with the infinite 
string coupling. The oj = —1 case is contained in this class. The low energy tree- level string 
action (a; = —1) with an axion leads to bouncing solutions as found in refs. [22[22l[22j, but 
these correspond to the solutions with divergent behavior of the dilaton and the curvature in 
past/future asymptotics. When the condition ()48|) is satisfied, the curvature asymptotically 
decreases but the string coupling diverges as r/ — > and ^ vr. Prom these arguments, we 
find that inclusion of the matter source (M ^ 0) is crucial for the construction of nonsingular 
bouncing solutions with finite dilaton and curvature even in the presence of the axion. In sec. 5 
we will analyze the dynamics of the system with and without axion in the context of string-gas 
cosmology. 

4.2 uj<-d/{d-l) 

Let us next proceed to the case of a; < —d/ (d—l). Introducing a new time parameter, rj, defined 

by 

* = ^,„(sinh,)^/<-"> with ^„=( -^'^"j;'-^ +"' )'"""', (49) 

the relation between r and rj yields 

dr 2 i+n 

:r = -7^00 (sin hry)!-" . (50) 

dr] (1 - n)\/E 
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In this case the value of N in eq. 1)411) can be either positive or negative^. 



4.2.1 Case of > 



The positive value of N can be realized in the presence of the modulus. Integrating eq. ()21() . 
the scale factor is given by the form ()25|) with 



tanh— 
2 




(51) 



Consider the evolution of a for q > with 77 > 0. When n < 1 both ^ and (p increase as r/ ^ 00, 
in which case the universe contracts asymptotically. In the case of n > 1, (j) decreases toward 
zero, whereas ^ approaches ^ — > ^0 as ry — > 00. In this asymptotic region, the scale factor evolves 
as 



a oc e("-i)(<«-i)'^ oc{t- to) with t 



(l + n)(d-l)-2d 

to + ce (i-^x-i-i) , 



(52) 



where to and c are constants. Since the power of the scale factor is larger than unity for n > 1, 
this represents an inflationary solution. Notice that this is achieved due to the presence of the 
modulus field. The curvature asymptotically decreases for 1 < n < whereas R diverges for 
n > In both cases the string coupling diverges in the limit of — > 00. The scale factor 

grows as the radiation dominant universe in another asymptotic, ^ — > -v/M. In this limit the 
scalar curvature decreases toward zero with a finite string coupling. Therefore it is possible to 
obtain nonsingular bouncing solutions for n > 1 if the universe starts out from the region with 
finite dilaton and curvature and evolves toward the weakly coupled regime. When the universe 
begins around the region ^ ~ ^/M and evolves toward </> = 0, we also have the bouncing solutions 
that asymptotically approach the inflationary solution H52() . However, since the solutions enter 
the strong coupling region in the latter case, this indicates the limitation of using the tree-level 
action. 

When q < the scale factor approaches zero as r/ — > 00 or 77 ^ 0, depending on the sign of 
q. This suggests a difficulty to obtain bouncing solutions for g < 0. In summary we can have 
regular bouncing solutions for q > and n > 1. 



4.2.2 Case of iV < 

When < the integration of eq. (PT|l leads to 



\N\ 



1 



M V^cos^ 



2/{d~l) 



(53) 



where 



^ = log Uo 



1 V 
tanh— 



In the absence of the modulus, A'' is negative. 




(54) 
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From this we find that bouncing solutions emerge in the interval rji < rj < rjf with 



Vi = log 



eo/^ + exp{(2j-l/2)7rM 



eo'/^-exp{(2j-l/2)^M 



Vf = log 



eo'/^ + exp{(2j + l/2)7rM 



exp{(2j + l/2)7r/p} 



Setting rj = rji + At] and r] 



ijf — Ar] around rj = rji,r]f, the scale factor evolves as 
a oc (A?7)"2/('^"^) (x{t- to)^ , 



(55) 



(56) 



where to is a constant. When n > 1 or n < 1/d, the solution represents the decelerating 
expansion of the universe with decreasing curvature {R 0). The standard axion coupling in 
d = 3 dimensions (n = —1) is included in this case, and the solutions exhibit the radiation 
dominant behavior with a oc (t — to)^^^- For the general axion coupling with n ^ —1, however, 
the evolution of the scale factor is not asymptotically radiation dominant. This means that the 
axion affects the dynamics of the system even in the asymptotic limit. In the case of l/d < n < 1, 
we have bouncing solutions with growing curvature in both asymptotics, which means that the 
solutions are not ideal ones. 

We have many sequences of intervals with r]i < rj < T]f, during each of which the string 
coupling and the curvature are kept finite. As in the case of < with no axion, nonsingular 
bouncing solutions can be constructed without entering the strong coupling region around ^ = 0. 



5 Application to string-gas cosmology 

In this section we shall consider the matter source term derived from the free energy of the ideal 
string-gas associated with type IIA/B string theory compactified on a square T^-torus, simple 
product of nine circles. It was recently shown that it is possible to construct the string-gas 
cosmological model where some dimensions (d) start to expand while the remaining dimensions 
are kept small ^7j- this work we include the effect of the axion. We also wish to investigate 
the validity of the analytic solutions in the presence of the large and small dimensions. In the 
region where the temperature, T = is below the Hagedorn one, Th = 2^/27r, the free energy 
of a string-gas in the canonical ensemble takes the following form |371 117j : 

where = a*^ is the volume of the large d dimensions with the scale factor a and 

A(r;f)= (58) 

represents the contributions of the whole Kaluza-Klein (KK) and winding modes along the small 
dimensions. Here f is the modular parameter of the world-sheet torus and q = exp(2i7rf). The 
sum in terms of p in eq. (|57|) runs only over positive odd numbers and corresponds to taking 
the correct quantum statistic for bosons and fermions. 

The M2 factor in eq. (|57j) is expanded in powers of g, as 

^2(r) = ^ = E D{Nm) , (59) 
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where 62 and f] are modular functions on the torus (see ref. pS!). Here the value of Nm corre- 
sponds to each string mass level with the degeneracy factor, D{Nm). When we integrate the ri 
and T2 integrals in eq. H57|) . it is convenient to extract the term with Nm = = rrii = rii = 
that corresponds to the purely massless state. Denoting this massless state (radiation) as -F^^f^, 
one has ^Jj 



Fial = -^m'r (^) (4^)(^+^)/^c(d+ i)(i - 2-('^-^i))r'^-^ , (60) 

where T and ^ are the Gamma function and the Riemann zeta-function, respectively. From this 
the energy E^'l^^ and pressure pj^^ along the d spatial dimensions are evaluated as 

^rad — J „^ ~ J ' v"^^ 





-df? 






= e^E, 




- V'A 


2 




= -e^P 
2 




-dfi^ 


2 




= -e^E 
2 



which represents the equation of state for radiation in d spatial dimensions. 

The pressure along the (9 — (i)-small dimensions is vanishing (i^^^^^ = 0) for the pure 
massless state. In this case the evolution of the small dimensions is trivial and can be kept 
small around the duality symmetric radius, r = 1 ^7]. Therefore it is sufficient to consider the 
evolution of only the large d-dimensions with a constant radius of the small dimension for the 
pure massless state. Introducing the shifted dilaton, ip = ip — dfj,, with a = e^, uj = —1 and 
X = 0, eqs. ©-© yield 

(63) 
(64) 

(65) 

Note that the energy and pressure are related with p and p in eqs. 0-®, as -E = 2a'^p and 
P = 2a'^p. Since the pure massless state satisfies the equation of state with p = dp from eq. 
the analytic method given in the previous section can be applied for this case. 

We are now considering the case of u; = — 1 and A = 0, thereby yielding a positive value of 
from eq. ()22() . In the absence of the axion we have q = ^^/3/6 from eq. ()24() . Since the negative 
value of q satisfies —1/2 < q < 0, bouncing solutions can be obtained in this case as shown 
in sec. 3.1 [see fig. This is achieved by choosing negative finite initial values of fi, in which 
case the scalar curvature is also finite. The scale factor exhibits a bounce at (p = (p^, given by 
(f^ . after which the system approaches the radiation dominant universe described by eq. (pTTj) 
as ^ — > (pc- In order to confirm these behaviors, we have numerically solved the background 
equations (jfv?|) -(|65 |) together with the equation for /3: 

Here we assumed an adiabatic evolution corresponding to the conservation of the entropy 5*. 
For the pure massless state (|6Up the inverse temperature, /3 = l/T, is proportional to the scale 
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factor o as in the standard radiation-dominated universe. We show in figs. [21 and |21 one example 
of the bouncing solution as dotted curves for the pure massless state without the axion . The 
dilaton continues to evolve toward the weakly coupled regime as expected, and the asymptotic 
behavior of the scale factor is found to be radiation dominant. 

Notice that the scale factor grows from the beginning for positive initial values of ji. In this 
case the time derivative of the dilaton, = ^/^ + d/i, can be positive even when ^ is negative 
initially. Since we are considering the weak coupling regime where e''' <C 1 is satisfied, it is typical 
that the r.h.s. of eq. ()63|) is negligible relative to the first and second terms. In the absence of 
the axion, ■0 is roughly estimated as ■0 ~ —y/dji. Then (p is positive for d > 1, thereby leading 
to the growth of the dilaton toward the strongly coupled regime. We argue that positive initial 
values of fi are not so welcome in order to avoid the breakdown of the string-gas description at 
finite temperature. In contrast, the bouncing solutions {i.e., /i < initially) correspond to the 
decreasing dilaton toward the weakly coupled regime. 

When the massive effect is taken into account, one has an extra source term in the small 
dimensions as well as the extra energy and pressure in eqs. (|6,3|) -(|65 j) . Let us study the leading 
terms that have an explicit dependence on the radius of the small dimensions (r = e'^) in the 
infinite sums appearing in eq. (|57|) . Taking the first KK and winding modes along a small 
direction with p = 1, i.e., the terms with {Nm = Nm = 0, mj = (1, 0, . . . , 0), nj = 0} (as well as 
rrii and ni exchanged), plus the remaining 8 — d inequivalent permutations, the pressure along 
the small dimensions is described as 



1 

^(d+3)/2-^("' 



(67) 



Here i^(d-i)/2 is the modified Bessel function, C{(3) is defined as C{I3) = (27r//3)('^+i)/2(18 - 
2c?)D(0)^/7r with -D(O) = 16 being a string degeneracy factor. We have an extra energy and 
pressure along the large dimensions due to the massive state in addition to (|61|) and (|62|) . See 
ref. ITj for the explicit forms of these terms. 

Numerically we solved the dilaton gravity equations (j63|) -(|65 |) and (j66|l together with the 
equation for the small dimensions with radius r = e^: 

v-^0 = ]^e^P^^-d) . (68) 

Since we are considering the case where the axion exists only for the large d dimensions, we do 
not include its effect along the small dimensions. The pressure completely vanishes at 

the self-dual radius r = 1 due to the compensation of KK and winding modes. Therefore the 
small dimensions are kept small around r = 1 provided that the initial values of v and v are 
close to zero. We have done numerical simulations for a variety of initial conditions and found 
that the massive effect for the large dimensions is weak as long as the small dimensions stay 
nearly constant around r = 1 (see figs.|21and|31as an example). Although the massive state gives 
rise to the extra energy and pressure, the dominant contribution comes from the pure massless 
state. For the ideal case where the small dimensions are kept small around r = 1, the system is 
well approximated by the state of pure radiation. 

Eliminating the energy E from eqs. ()63|) and H65|l and making use of the approximate relation 
E ~ dP in eq. (|64)) . we have 

1 1 1 p(n+l)v 
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Figure 2: The evolution of the scale factor a in the case of d = 3 with or without the axion for 
the initial conditions fiQ = —0.03, fio = 3.0, z>o = —0.005, fo = 0.01, ipo = —14, Pq = 17. The 
massive effect is taken into account together with the radiation for the solid curves, whereas the 
dotted curve corresponds to the pure radiation case without the axion. 
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Figure 3: The evolution of the dilaton 99 in the case of d = 3 for the same initial conditions 
shown in fig. |2j The massive effect is taken into account together with the radiation for the solid 
curves, whereas the dotted curve corresponds to the pure radiation case without the axion. 
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Figure 4: Phase space trajectories of the solutions of string-gas cosmology for d = 3. We 
implement the effect of the massive state together with the leading massless state. In the 
absence of the axion (solid curves), the solutions asymptotically approach the line given by 
l/f = — which represents the radiation dominant universe. The expansion rate I gets larger 
when the axion is included (dotted curves), but the solutions eventually approach the radiation 
dominant universe. 



2d 2d a?- 



i - -7;l' + fl + iTJ+ ' :J B^-^^ (70) 



where / and / are defined as f = ip and X = I. In the absence of the axion (B = 0), it was shown 
in refs. jl9l I17j that the line / = —f/d is an attractor for the massless state. From fig.|l]we find 
that the solutions without the axion (solid curves) asymptotically approach the line I = — f/d 
even in the presence of the massive state. 

When the axion is included in the large d dimensions, the evolution of the scale factor and 
the dilaton should be changed as discussed in sec. 3.1. As one example we show in figs. HI and |31 
the evolution of a and (p for d = 3 and n = —1. The initial conditions are chosen to he ipi < 
and e"^* <C 1. The dilaton does not continue to evolve toward the weakly coupled region [see 
eq. ()42() ] except for the initial stage, unlike the case without the axion. In fact, we find from fig-IHl 
that the dilaton changes to evolve toward the strong coupling regime after the initial decrease. 
The presence of the axion provides additional source terms in eqs. (|69() and ()7U() . thereby leading 
to the different evolution of the scale factor as shown in fig. |2j When n = —1, the growth of 
the scale factor is faster than in the case without the axion due to the positive last term in the 
r.h.s. of eq. (|70)) . Nevertheless we find from fig. |3 that the solutions eventually approach the 
line I = — f/d as / and / decrease toward zero. This means that the dilaton approaches a finite 
value with the radiation dominant universe (a oc as ^ — > = ^f^), as discussed in 

sec. 4.1. In addition the scalar curvature decreases to zero as ^ — > ^c- 

Since we are now dealing with the case of w = — 1 and ^ = 0, the condition (|46() is satisfied 
whereas (|48j) is not. The condition (|48|) corresponds to the bouncing solutions where the scalar 
curvature R converges in the limit of (/> — > 0, whereas R diverges in the present case in this limit 
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{<jj = —1,^ = 0). Nevertheless, when the universe starts out from a weak couphng regime with 
a finite curvature, it is possible to construct bouncing solutions that asymptotically approach 
the radiation dominant universe with finite R and (/? irrespective of the presence of the axion 
as found in sec. 4.1 for the pure massless case. The examples plotted in figs. Hand 01 show the 
existence of nonsingular bouncing solutions even when the condition H48|) is not fulfilled. Notice 
that in the presence of the axion but without the radiation-like source matter ||21j R and (p 
diverge in both asymptotic limits. This indicates that inclusion of the matter source such as the 
ideal string-gas is crucial to construct the bouncing solutions without singularities in asymptotic 
future. It will be interesting to extend our analysis to other matter sources which appear e.g., 
in the p-brane (Dp-brane) cosmology. 

We also made numerical analysis for larger initial values of the dilaton which are close to 
~ 1. We found that the system tends to evolve toward the strongly coupled regime with 
the growth of toward infinity. The scalar curvature shows divergent behavior together with 
the unbounded growth of the dilaton. Nevertheless we have to keep in mind that the string- 
gas description based on the free energy (|57]) ceases to be valid when the system enters the 
regime of strong coupling and large curvatures. In addition we need to consider the higher-order 
corrections to the tree-level effective action in this regime. In the next section we will discuss 
the dynamics of the system when the dilatonic higher-order corrections are taken into account. 



6 Inclusion of higher-order corrections to the tree-level action 

In order to understand the effect of the higher order corrections to the tree-level action, let us 
consider the dilaton-dependent corrections, given by (Sj 

Cc = -\c^'K{v) [ciRIb + C2(Vv9)^] , (71) 

where 5,{ip) is a general function of (p. A is an additional parameter depending on the species of 
strings. For example, we have A = —1/4,-1/8 for the Bosonic and Heterotic string, respectively. 
The Gauss-Bonnet term, iJ^g = R^ — AR^'" R^v + R^'^^^Rnua/s, has the property of keeping the 
order of the gravitational equations of motion unchanged. The dilatonic corrections Cc are the 
sum of the tree-level a' corrections (i = 0) and the quantum i-loop corrections (i = 1, 2, 3, • • •), 
with the function ^((/j) given by [51 [TUl 1^ 

^(<^) = - ^ ae(-i)^ , (72) 

n=0 

where Ci (i > 1) are the coefficients of i-loop corrections with Co = 1. 

If the universe starts out from a sufficiently weak coupling regime {e'^ <C 1), it can reach 
the high curvature regime while the string coupling is still small. In this case the a' correction 
(n = 0) becomes important before the quantum loop corrections begin to work. Gasperini et 
al. jS| showed that inclusion of the a' correction gives rise to the fixed points in the {(p, H) 
plane when only the dilaton is present. This leads to an inflationary solution driven by the 
constant H and the growing (p. Although the successful transition from this dilaton-driven 
inflationary phase to a decelerated FRW era is difficult to be achieved only in the presence of 
the a' correction, this problem can be solved by implementing the quantum loop corrections. In 
fact Brustein and Madden [Hj showed that inclusion of the quantum loop corrections can lead 
to the successful graceful exit due to the violation of the null energy condition. 
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Here we wish to analyze the dynamics of the system in the presence of the axion, modulus 
and the stringy matter source. Let us consider the case of d = 3 and uj = —1. The source terms 
due to £c in eqs. 0-® are given bv P^ITU] 



where 



{Pc}i 
{Pch 

{Ac}, 



i=0 i=0 i=0 



a 



(73) 



a 



I'X^e^'-^^'^ {-24ci(i - l)H^ + 3c29i^} , (74) 

y'\ei^-^)'P |8ci(i - l)H [(i - l)^'^H + (pH + 2^{H + H^)\ + cz^i^} , (75) 

/Ae(^-^)^ |24ci(i - l)H^{H + F^) - 3c20^ [4^3 + 4vi/? + (i - 1)93^] | . (75) 

Substituting these expressions for eqs. ©-©) we get the equations for H and ip : 



6 - 24a'AciF2 ^ C7,(i - l)e*'^ 



j=0 



(77) 



-ne 



(n-i)^^2 ^ ^/^ ^ (^.gi<p |24ci(i - l)i/^ - 3c2vi2[4(^ + 4Hip + {i - l)p^]} 



i=0 



2-8a'Aci-H'2^Ci(z- l)e*'^ 



i=0 



■ 2 



(78) 



+e-(n-iV^2 _ a'Xj2Qe'^ {8ci(i - l)H[{i - l)^'^H + 2ip{H + i/^)] + C20^} + 2e^p . 

In the case of x = o" = and p = p = 0, the solutions approach the string phase characterized by 
constant H and ip provided that the system begins from the weakly coupled regime {e'^'- <C 1). 
These values are obtained by setting H = and (p = in eqs. ((77|) and ((75|) with i = 0. The 
effect of the quantum loop corrections (i 7^ 0) begin to work around e'^ ^ 1, which can lead to 
the graceful exit toward the declerating FRW branch. The number of e-folds, iVg = In a, before 
the graceful exit is approximately estimated as ~ \ipi\/2, where \(pi\ is the initial value of the 
dilaton. 

Let us study the case ofa; = — 1,X7^0 ^-^d cr / 0, and p = p = 0. In the dilaton-driven 
phase corresponding to the lowest order PBB cosmology, the evolution of a and (j) is described 
by a oc (— i)^/^ and (/) = oc (— t)^"*"^. From eq. (fTOj) x and a evolve ds, x cc {—t)~^ and 
a oc (— t)v^~"(v^+i). Therefore x grows during the dilaton-driven phase, whereas a decreases 
for n < (3 - ^/3)/2 . The growth of X can change the dynamics of the system through eqs. (|77j) 
and (|78|) . For a fixed value of the initial Hubble rate, the initial 'p gets larger from the constraint 
equation in the presence of the modulus. We also find from eq. (|78|) that the increase of x 
leads to the faster growth of (p. The axion has a similar effect for n > (3 — \/3)/2. 

The growth of x (and a for n > (3 — \/3) /2) terminates once the system enters the string 
phase characterized by the constant H and <p {Hq = 0.616 and ipo = 1.40 for d = 3). During this 
phase, we have a oc e^"* and (p oc (pot, thereby yielding x oc e-i3^o-'Po)i g^^id a oc g-CSJ^o-'T.v'o)*^ 
The fixed point {Hq,ipq) exists in the range of 3Hq > ipQ in order to allow for a subsequent 
graceful exit j40j . Therefore x begins to decrease once the system enters the string phase, 
and the situation is the same for a as long as n < SHo/ipQ. In standard low-energy effective 
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Figure 5: The evolution of H and 99 as a function of the number of e- folds for d = 3, p = p = 
and n = —1 with initial conditions ipi = —130 and Hi = 1.5 x 10~^. We include the quantum 
corrections up to i = 2 with coefficients Ci = 1 and C2 = — 10~^ together with the a' correction 
(Co = 1). The dotted curves correspond to the evolution without the modulus and axions 



fields, whereas these fields are involved in the solid curves with initial conditions Xi 
di = 10~^. We do not find a significant difference even when x and a are present. 
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Figure 6: The evolution of the dilaton (p in the presence of the modulus and axions fields for 
the same initial conditions as in fig. [21 The initial dilaton-driven phase (0 < t < 100) is followed 
by the string phase with a linearly growing dilaton (100 < t < 200), after which the solutions 
connect to the declerating FRW branch. Inset: The evolution of x Both x a are 

exponentially suppressed during the string phase. 
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string theory we have n = — 1 for the axion couphng, in which case a is exponentially 
suppressed during the string phase and its energy density becomes very small relative to the 
dilaton. Typically n is smaller than unity not only in low energy effective string theory but 
also in multidimensional models with gauge fields or in the higher dimensional conformal gauge 
field theory j24j . In such cases the effect of the axion is not dynamically important when the 
correction (f7T|) is present. 

We numerically solved the equations (|77|) and (|78|) for uj = —\ and n = —1 together with the 
initial conditions satisfying e"^"- ^ 1, < 1 and o"j < 1. In the presence of the modulus and the 
axion, the string phase with constant H and ip is typically reached earlier than in the case only 
with the dilaton. Nevertheless the difference is not significant as seen from fig. |S1 In both cases 
the quantum loop corrections begin to work around e*^ ~ 1, after which the solutions connect 
to the declerating FRW branch. When only the a' and one-loop corrections are present (Co = 1 
and Ci = 1), the dilaton and the curvature exhibit unbounded growth after the branch change 
as shown in ref. j9j in the single field case. If the two-loop correction with negative coefficient 
C2 is included in eq. (|72|1 , this can overwhelm the one-loop correction when the dilaton becomes 
sufficiently large. In fact, as shown in fig. El the evolution of the dilaton becomes mild by 
including the two-loop correction. 

Even when the matter source is present {p 7^ 0,p 7^ 0), the dynamics of the system is 
dominated by the correction H71|) . The effect of the matter is involved in the last term of 
eq. (|78|). Since we are considering the initial conditions corresponding to the weakly coupled 
region (e'^^ <^ 1), the a' correction tends to dominate over the last term in eq. H78|) . Then the 
system soon approaches the string phase with constant H and (p. During the string phase, p 
and p are rapidly suppressed due to the exponential growth of the scale factor. This means 
that the matter source term is practically negligible around the graceful exit. For example, the 
pressure of the string-gas discussed in the previous section evolves as p oc oc e"^*^ in the 
string phase, thereby leading to the exponential decrease of the last term of eq. (f78|) . In fact we 
have done numerical simulations for the free energy (|6U() in the context of string-gas cosmology 
and found that the evolution of the system is almost the same as shown in figs. El and El The 
energy density and pressure of the string-gas are exponentially suppressed, together with the 
rapid decrease of the temperature (T oc a~^) before the graceful exit. 

We also investigated the case where the dilatonic Brans-Dicke parameter to does not equal 
to —1. When uj lies in the range —3/2 < uj < with d = 3 and A = 0, the qualitative behavior 
is mostly the same as in the a; = — 1 case. The difference is that the values of H and (p in 
the attractor region get bigger for smaller \io\, thereby leading to larger amount of e- foldings. 
However, when lo < —3/2, we find that the system is typically unstable and tends to fall into 
a singularity. Although nonsingular solutions can be obtained when uj is close to —3/2 and e"^ 
is not too much smaller than unity, the parameter range of these solutions is restricted to be 
very narrow. This suggests that inclusion of the correction Cc is not necessarily welcome for the 
construction of regular solutions for a; < —3/2. 

7 Summary and discussions 

In this paper we have studied the construction of nonsingular cosmological solutions in a gen- 
eral Z)-dimensional effective action derived from string theories. In particular we tried to find 
singularity-free solutions analytically for the dilaton-modulus-axion system in the presence of 
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the "stringy" matter source term. We also keep the dilatonic Brans-Dicke parameter u arbitrary 
so that the action contains a variety of theories such as the F-theory and the multidimensional 
theory. In the tree- level action without higher-order a' and loop corrections, we found the 
following results. 

• Even for the general {d+ l)-dimensional action with modulus and axion fields, the system 
is integrable for the flat FRW background if the matter source is radiation (p = dp) or 
there is no matter {p = p = 0). 

• In the case of = —1, we find nonsingular bouncing solutions where the universe starts 
out from the state of the finite scalar curvature and asymptotically approaches the radia- 
tion dominant stage with decreasing curvature and finite string coupling. This situation 
naturally appears in the context of string-gas cosmology. 

• In the theories where a; is largely negative, we have a sequence of regular bouncing solutions 
where the string coupling and the scalar curvature remain finite. This was already pointed 
out in ref. [2^1 EEl for (i = 3, but we extended the analysis to the case of the general spatial 
dimension d with modulus fields. 

• As a "stringy" matter source, we considered the ideal string-gas in thermal equilibrium 
taking into account the massive state coming from Kaluza-Klein and winding modes in 
addition to the leading massless state. Numerically we obtained nonsingular bouncing 
solutions where the "large" 3 dimensions asymptotically approach the almost radiation 
dominant universe irrespective of the presence of the axion. This indicates that inclusion 
of the radiation-like source is crucial for the existence of such regular solutions. 

While bouncing solutions can be obtained in the tree-level action with the radiation-like 
matter source, it can happen that the system enters the strongly coupled region. In this case 
it is inevitable to implement the higher-order corrections to the tree-level action. Therefore 
we also analyzed the case where the dilatonic higher-derivative and loop corrections are taken 
into account for the dilaton-modulus-axion system with some matter source. We find that 
the dynamics of the system for w = — 1 is not significantly changed by including scalar fields 
or matters other than the dilaton. The system approaches the string phase characterized by 
constant H and 99, after which the graceful exit is realized through the dominance of higher-order 
loop corrections. Typically the energy densities of the modulus, axion, and the matter source 
are exponentially suppressed once the dilaton-driven infiation begins in the string phase. We 
also find that the system tends to be unstable for the case of large negative oj. 

There are several open issues that we do not study but that deserve further investigation. 
While we consider the fiat FRW background for simplicity, it should be fairly easy to extend 
our analysis to the closed/open FRW background with general D = d + 1 dimensions, as was 
already done for the case of d = 3 without the modulus ^O]- In addition it is worth analyzing 
the anisotropic background in order to check the generality of the singularity avoidance. In fact 
a new singularity- called the determinant singularity- appears in the presence of anisotropics 
when the action contains a Gauss-Bonnet term jlj. It is of interest to investigate whether the 
singularity solutions we found are stable by including small anisotropies. 

The validity of our nonsingular cosmological scenarios may be investigated further by an- 
alyzing the power spectra of produced density perturbations. Recent observations of WMAP 



24 



j42j suggest that the spectrum of scalar perturbations is close to scale-invariant. In the sim- 
plest tree-level pre-big-bang scenario, the curvature perturbation TZ is blue-tilted with a spectral 
index n ~ 4 43 , which is incompatible with observations. This situation does not change so 
much even involving the higher-order corrections to the tree- level action |44( I45j . Nevertheless 
it is possible to have nearly scale-invariant spectra of the axion perturbations depending on the 
expansion rate of the modulus field j22j. Making use of the "curvaton" mechanism j42j, there 
remains a possibility to explain the observationally supported flat spectra if the axion plays the 
role of the curvaton (see ref. In addition, originating from the proposal of Ekpyrotic cos- 

mology there are a number of controversial argument j48| llOj about the choice of the gauge 
conditions for scalar perturbations generated in bouncing cosmologies. It is certainly of interest 
to analyze the evolution and the spectra of scalar perturbations in our nonsingular bouncing 
models by taking into account the curvaton perturbations. 
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